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Abstract. We present the results of the first search for gravitational wave bursts associated 
with high energy neutrinos. Together, these messengers could reveal new, hidden sources 
that are not observed by conventional photon astronomy, particularly at high energy. Our 
search uses neutrinos detected by the underwater neutrino telescope ANTARES in its 5 line 
configuration during the period January - September 2007, which coincided with the fifth 
and first science runs of LIGO and Virgo, respectively. The LIGO-Virgo data were analysed 
for candidate gravitational-wave signals coincident in time and direction with the neutrino 
events. No significant coincident events were observed. We place limits on the density of 
joint high energy neutrino - gravitational wave emission events in the local universe, and 
compare them with densities of merger and core-collapse events. 
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1 Introduction 


Multi-messenger astronomy is entering a stimulating period with the recent development of 
experimental techniques that will open new windows of cosmic radiation observation in all its 
components. In particular, both high-energy (3> GeV) neutrinos (HENs) and gravitational 
waves (GWs), neither of which have yet been directly observed from astrophysical sources, 
are becoming new tools for exploring the Universe. 

While HENs are expected to be produced in interactions between relativistic protons 
and the external radiation field of the source (e.g., [85, 146]), GWs carry information on the 
intricate multi-dimensional dynamics in the source’s central regions (e.g., [58]). HENs and 
GWs are thus complementary messengers. 

Simultaneous emission of GWs and HENs has been proposed in a range of cataclysmic 
cosmic events including gamma-ray bursts (GRBs), core-collapse supernovae (CCSNe), soft- 
gamma repeater outbursts (SGRs), and, potentially, cosmic string cusps in the early universe. 

Observational constraints on HEN and GW emission from some of these phenomena 
have already been obtained. The IceCube collaboration recently placed limits on the HEN 
emission in GRBs [8, 10, 11], SGRs and blazars [13], and jet-driven CCSNe [12] using data 
from the IceCube detector at various levels of completion. Similarly, the ANTARES Collab- 
oration has placed limits on the HEN flux from gamma-ray flaring blazars [25] and GRBs [24], 
as well as a diffuse muon neutrino flux from extragalactic sources [33]. These limits, however, 
do not yet strongly constrain HEN emission and ultra-high-energy cosmic-ray production in 
relativistic outflows [80, 85, 99]. The LIGO Scientific Collaboration and Virgo placed limits 
on the GW emission in GRBs [1, 6, 20] and SGRs [4, 15, 19]. The exclusion distances of 
these searches were, however, not sufficiently large to expect a GW detection. 

The above HEN and GW searches used timing and sky location information from ob- 
servations of events in the electromagnetic spectrum. A potentially large subset of GW and 
HEN sources may be intrinsically electromagnetically faint, dust-obscured, or missed by tele- 
scopes, but sufficiently luminous in GWs and HEN to be detected. Such sources may include, 
but are not limited to, partially or completely choked GRBs with, perhaps, only mildly rel- 
ativistic jets [38, 113, 130, 143], relativistic shock breakout in compact CCSN progenitor 
stars [147], and cosmic string cusps [45, 59, 134]. 

Searches for HENs and GWs from such events have thus far relied on blind (i.e., un- 
triggered) all-sky searches. An all-sky search for point sources of HENs in IceCube data was 
performed [13] and a similar study was carried out with ANTARES data [26]. LIGO and 
Virgo have carried out a number of all-sky searches for GWs. The most recent and most sen- 
sitive such search for model-independent GW bursts was published in [5], whereas the most 
recent allsky search for binary inspiral-merger can be found in [7]. All-sky model-dependent 
constraints on cosmic string GW emission have been placed [17]. The sensitivity of such 
blind all-sky searches is limited by a much larger background compared to searches based on 
timing and sky locations from electromagnetic observations. 

A search for temporally and spatially coincident HEN and GW signals is a strong al- 
ternative to electromagnetically triggered or blind all-sky analyses that search for GWs or 
HENs individually. Such a search is independent of bias from electromagnetic observations, 
but still enjoys a much reduced background thanks to timing and sky location constraints. 
A similar idea has been used in the follow-up of candidate GW events by the low-energy 
neutrino detector LVD [29]. A joint GW-HEN search was first proposed in [39] and [127], 
and constraints on joint GW-HEN signals based on the interpretation of independent GW 
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and HEN observational results were derived in [43]. Here we present the first direct search 
for coincident GW-HEN events, using data taken by the ANTARES HEN telescope and 
by the LIGO and Virgo GW observatories from January to September 2007. At this time, 
ANTARES was still under construction and operating with only 5 active lines. At the same 
time, the fifth LIGO science run (S5) and the first Virgo science run (VSR1) were carried 
out. This was the first joint run of the LIGO-Virgo network with the detectors operating 
near their design sensitivities. 

The basic principle of the analysis presented here is that of a “triggered” search: HEN 
candidates are identified in the ANTARES data, then the GW data around the time of the 
HEN event are analyzed for a GW incident from the HEN estimated arrival direction. This 
method has been applied previously in searches for GWs associated with GRB triggers [6, 20]. 
It has been shown to have a distance reach up to a factor of 2 larger [6] than a blind all-sky 
search of the GW data, due to the reduced background. The expected rate of detections 
depends also on the beaming of the trigger signal, since the triggered search is only sensitive 
to the subset of sources oriented towards Earth. The comparison of the analysis method 
used in this paper to the LIGO-Virgo blind all-sky search [5] has been done in [144], and 
predicts a detection rate for the triggered search of between 0.1 and 6 times that of the 
blind search for beaming angles in the range 5°-30°. These numbers are broadly consistent 
with estimates for the special case of dedicated matched-filter searches for compact binary 
coalescence signals associated with short GRBs [e.g., 54, 63, 89, 116] after rescaling for a 
smaller distance improvement factor (typically ~1.3, due to the better inherent background 
rejection of these specialised searches). In either case, most of the GW events found by the 
triggered search will be new detections not found by the all-sky blind search, illustrating the 
value of the triggered search even when the relative detection rate is low [144], 

We analyze a total of 158 HEN events detected at times when two or more of the LIGO- 
Virgo detectors were operating. ANTARES is sensitive to HENs with energies greater than 
~ 100 GeV [27]. The LIGO-Virgo analysis targets model-independent burst GW signals with 
durations < 1 s and frequencies in the 60 Hz to 500 Hz band. The GW search is extended in 
frequency up to 2000 Hz only for a subset of the HEN events, because the computational cost 
of such a search with the current GW analysis pipeline is prohibitive. 

Statistical analyses of the HEN sample show no sign of associated GW bursts. 

This paper is organized as follows. In section 2 we discuss sources of coincident HEN 
and GW emission and expected prospects for their detection. In section 3 we describe the 
ANTARES, LIGO, and Virgo detectors and the joint data taking period. Section 4 de- 
scribes how the HEN sample was selected. Section 5 describes the search for GWs coincident 
in time and direction with the HEN events. We present the results of the search in sec- 
tion 6. We discuss the astrophysical implications of the results in section 7 and conclude 
with considerations of the potential for future joint GW-HEN searches. 

2 Candidate sources for high-energy neutrino and gravitational wave emis- 
sion 

HEN emission is expected from baryon- loaded relativistic astrophysical outflows. In the most 
common scenario (e.g., [146]), Fermi-accelerated relativistic protons interact with high-energy 
outflow photons in p'j reactions leading to pions or kaons, whose decay results in neutrinos, 
e.g., 7T + — > p, + + Vp — > e + + v e + v f j-\-v il , which is the dominant process (see, e.g., [149]). This 
gives {y e + v e : i/^ + P /t : v T + u T ) production ratios of (1 : 2 : 0), changing to approximately 
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(1 : 1 : 1) at Earth due to flavor oscillations (e.g., [40]). The HEN spectrum depends on 
the spectrum of the accelerated protons (e.g., [8, 76, 85]) and, thus, on the properties of 
the astrophysical source. In this section, we provide estimates of the sensitivity of the 5- 
line ANTARES detector for HENs from the various potential sources by estimating the 
probability P — X% that at least one HEN is detected for a source at a given distance d\- 

GW emission occurs, at lowest and generally dominant order, via accelerated quadrupo- 
lar mass-energy dynamics. The coupling constant in the standard quadrupole formula for 
GW emission (e.g., [140]) is Gc~ 4 ~ 10" 49 s 2 g" 1 cm -1 and the directly detectable GW strain 
scales with (distance)" 1 . For example, a source at lOMpc needs a quadrupole moment of 
1M 0 x (100 km) 2 that is changing on a millisecond timescale to be detectable by a GW 
detector sensitive to a strain of 10“ 21 . Equivalently, the minimum GW energy emission de- 
tectable by the LIGO-Virgo network at this distance is approximately -Egw — 10 “ 2 M 0 c 2 
to 10 MqC 2 for frequencies between 100 Hz and 1000 Hz [5]. 

In the following, we discuss a number of astrophysical scenarios in which both HENs 
and GWs may be emitted at detectable levels. 

2.1 Canonical long gamma-ray bursts from massive stars 

Long-duration GRBs (LGRBs; Tgo > 2 s; Xgo is the time over which 90% of the 7 counts 
are detected) are observationally implicated to be related to the collapse of massive stars 
normally leading to core-collapse supernova explosions [83, 110]. Typical LGRBs are strongly 
beamed and most likely have jets with Lorentz factors T > 100 and isotropic equivalent 
luminosities of 10 51 ergs" 4 to 10 53 ergs _1 [68, 107, 123]. LGRBs are detected at a rate of 
order l/(few days) by 7-ray monitors on satellite observatories such as Swift /BAT [42, 69] and 
Fermi/GBM [46, 105]. It is important to note, however, that Swift/BAT and Fermi/GBM 
miss ^90% and ^40% of the prompt emission of all GRBs, respectively. This is due to limited 
fields of view, technical downtime, and orbital passes through the South Atlantic Anomaly. 
A nearby (say ^lOMpc) LGRB will have a bright multi-wavelength afterglow and may be 
accompanied by a CCSN (see section 2.2), but a significant fraction of local CCSNe may 
have been missed by CCSN surveys based on galaxy catalogs [104]. 

The central engine of LGRBs is expected to either be a collapsar (a black hole with 
an accretion disk; [103, 150]) or a millisecond magnetar (an extremely rapidly spinning, 
extremely magnetized neutron star; e.g., [109]). In both scenarios, HEN emission may result 
from a relativistic expanding fireball. HENs may begin to be produced even before the jet 
breaks out of the stellar envelope [128] and may continue well into the afterglow phase [112]. 

HEN emission from canonical LGRBs is expected to have appreciable flux for energies 
in the range 100 GeV to 100 TeV. For a LGRB at ^50 Mpc (^lOMpc) one would expect 
of order 1 (100) HEN events in a km 3 -scale water- or ice-Cherenkov detector (e.g., [8, 76, 
85, 146]). Based on the flux predictions of [76], the probability for detection in the 5-line 
ANTARES detector can be estimated to be ^50% for a source at 10 Mpc, which decreases 
to ^2% for a source at 50 Mpc. Note that these are most likely optimistic estimates, since 
more detailed analyses suggest lower HEN fluxes from GRBs [e.g., 85]. 

The most extreme scenario for GW emission in LGRBs is dynamical fragmentation of 
a collapsing extremely rapidly differentially spinning stellar core into a coalescing system of 
two protoneutron stars [61, 93]. Such a scenario may be unlikely given model predictions for 
the rotational configuration of GRB progenitor stars [e.g., 151]. Its GW emission, however, 
would be very strong, leading to emitted energies Fqw ~ 10~ 2 M 0 c 2 to 10" 1 M 0 c 2 in the 
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50 Hz to 1000 Hz frequency band of highest sensitivity of the initial LIGO/Virgo detectors, 
which could observe such an event out to approximately 20Mpc to 40Mpc [6, 7]. 

In more moderate scenarios backed by computational models, GW emission from LGRBs 
is likely to proceed, at least initially, in a very similar fashion as in a rapidly spinning 
CCSN [67, 95, 118]. If a black hole with an accretion disk forms, a second phase of GW 
emission may come from various hydrodynamic instabilities in the accretion disk [e.g., 92, 
124, 141], 

In the initial collapse of the progenitor star’s core, a rapidly rotating protoneutron star 
is formed. In this process, a linearly-polarized GW signal is emitted with typical GW strains 
\h\ ~ ICG 21 to 10 -20 at a source distance of lOkpc and emitted energies £gw ~ 10 “ 8 M 0 c 2 
to 10“ ' Mqc 2 between 100 Hz and 1000 Hz [64, 120]. This part of the GW signal will only 
be detectable for Galactic events and is thus not relevant here. 

In its early evolution, the protoneutron star (or protomagnetar, depending on its mag- 
netic held) may be spinning near breakup. This can induce various rotational instabilities that 
induce ellipsoidal deformations of the protoneutron star, leading to strong, quasi-periodic, 
elliptically-polarized GW emission [55, 66, 67, 118, 131]. A typical GW strain for a deformed 
protoneutron star of 1.4 M 0 and radius of 12 km, spinning with a period of 1ms may be 
h ~ few x 10~ 22 at lOMpc. If the deformation lasted for 100 ms, £gw ~ 10 -1 M 0 c 2 would 
be emitted at 2000 Hz [66]. 

In the collapsar scenario, accretion onto the protoneutron star eventually leads to its 
collapse to a spinning black hole [117]. This and the subsequent ringdown of the newborn 
black hole leads to a GW burst at few x 10 2 Hz to few x 10 3 Hz with h ~ 10 -20 at 10 kpc 
and £gw ~ 10“ 7 M q c 2 . It is thus detectable only for a Galactic source [119]. 

More interesting are hydrodynamic instabilities in the accretion disk/torus that forms 
after seconds of hyperaccretion onto the newborn black hole. The inner parts of the disk are 
hot, efficiently neutrino cooled and thus thin [e.g., 126] while the outer regions are inefficiently 
cooled and form a thick accretion torus. Gravitational instability may lead to fragmentation 
of this torus into one or multiple overdense regions that may could condense to neutron-star- 
like objects and then inspiral into the central black hole [124]. For a source at lOMpc, a 
IMq fragment and a 8 M 0 central black hole, this would yield strains of h ~ few x ICG 22 and 
emitted energies in the most sensitive band of 10 3 Mqc 2 to 10 2 Mqc 2 . 

The accretion torus may be unstable to the Papaloizou-Pringle instability or to co- 
rotation-type instabilities [121, 122], h ~ 10~ 21 to 10 -2 ° was estimated in [92] for a source at 
10 Mpc and GW frequencies of 100 Hz to 200 Hz for am — l — dominated non-axisymmetric 
disk instability in a disk around a 10 Mq black hole. This corresponds to -Egw of order 
10 -2 M 0 c 2 to 10 _1 Mqc 2 . 

In the speculative suspended accretion model for GRB accretion disks [141], low-order 
turbulence powered by black-hole spindown may emit strong GWs. In the frequency domain, 
this results in an anti-chirp behavior, since most of the emission is expected to occur near 
the innermost stable orbit, which moves out in radius as the black hole is spun down. Simple 
estimates suggest GW strains h ~ HR 21 at 10 Mpc and frequencies in the 100 Hz to 1000 Hz 
band. Depending on the initial black hole spin, £gw could be of order 1 M 0 c 2 . 

2.2 Low-luminosity GRBs and engine-driven supernovae 

Low-luminosity GRBs (LL-GRBs; also frequently referred to as X-ray flashes) form a subclass 
of long GRBs with low 7-ray flux (e.g., [57, 83, 110]). LL-GRBs are much more easily 
missed by observations than LGRBs (see section 2.1) and the small observable volume (due 
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to their low luminosity) suggests an event rate that may be significantly higher than the 
rate of canonical LGRBs [52, 97, 100, 135, 142], Five of the seven GRBs that have been 
unambigously associated with a CCSN are LL-GRBs [83, 110, 152]. Moreover, all GRB- 
CCSNe are highly energetic and of the spectroscopic type Ic-bl subclass. Ic indicates a 
compact hydrogen/helium poor progenitor star and the postfix -bl stands for “broad line,” 
because they have relativistically Doppler-broadened spectral features. 

Type Ic-bl CCSNe occur also without LL-GRB or LGRB, but are frequently identified 
as engine-driven CCSNe that exhibit luminous radio emission [e.g., 135, 136]. 

Theory suggests (e.g., [50, 51, 90, 96]) that the transition between engine-driven CCSNe, 
LL-GRBs with CCSNe, and canonical LGRBs may be continuous. All are likely to be driven 
by a central engine that launches a collimated bipolar jet-like outflow and their variety may 
simply depend on the power output and duration of central engine operation [50, 96]. The 
power output of the engine determines the energy of the jet and its Lorentz factor. The 
duration of the central engine’s operation determines if the jet can leave the progenitor star 
and make a normal LGRB. If it fails to emerge, the LGRB is “choked” and a more isotropic 
energetic CCSN explosion is likely to result. As suggested by [50], the relativistic shock 
breakout through the stellar surface could then be responsible for a LL-GRB. 

The GW emission processes that may be active in LL-GRBs and engine-driven CC- 
SNe are most likley very similar to the LGRB case discussed in section 2.1 and we shall 
not consider them further here. HEN emission is expected from the entire range of stellar 
collapse outcomes involving relativistic flows. Since LL-GRBs and engine-driven CCSNe are 
most likely much more frequent that canonical LGRBs, much effort has been devoted to 
understanding the HEN emission from such events [38, 84, 108, 112, 113, 128-130, 143, 147]. 
It is worthwhile to consider the probability of detection of HENs in the 5-line ANTARES 
detector from LL-GRBs and engine-driven CCSNe, in which mildly relativistic jets are likely 
to be involved. The detection probability depends strongly on the energy and the Lorentz 
factor of the jet. Using the reference parameters of [38], T = 3 and £j e t = Lo ~ 3 x 10 51 erg, 
the detection probability is ^50% at dso = 1 Mpc. 

2.3 Mergers and short gamma-ray bursts 

Short-duration GRBs (SGRBs; Tqq < 2 s) are rarer than LGRBs and expected to result 
from double neutron star (NS-NS) and/or neutron star - black hole (BH-NS) mergers [e.g., 
68, 114]. The isotropic equivalent energy of SGRBs is 2 to 3 orders of magnitude smaller 
than the energy of LGRBs. Their jets have most likely lower Lorentz factors of T ~ 10 to 50 
and wider opening angles. Due to their short duration and low isotropic equivalent energy, 
SGRBs are much easier to miss observationally than LGRBs and their observable volume is 
much smaller. 

The efficiency of HEN emission in SGRBs depends on the efficiency of proton accelera- 
tion, the 7-ray flux, and the SGRB variability time scale [114]. For a simple estimate of the 
detection probability in the 5-line ANTARES detector, one may resort to the HEN flux es- 
timates of [76] (but see [85] for refined results). Assuming a jet with T = 300, E- ]e t = 2x 10 50 
erg, one finds a distance dio ~ 10 Mpc at which the probability of HEN detection by the 5- 
line ANTARES detector is 10%. Hence, only the closest and/or the most powerful SGRBs 
may be detectable. 

The GW emission from NS-NS and BH-NS mergers is well studied [see 65, 133, for 
reviews]. Most of the emission comes from the late inspiral and merger phase during which 
the binary sweeps through the 50 Hz to 1000 Hz band of highest sensitivity of LIGO/Virgo. 
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The total emitted -Egw is of order 10“ 2 Mqc 2 to 10 _1 Mqc 2 . At the time of this analysis the 
LIGO/Virgo network had maximum sensitive distances of ^30Mpc for equal-mass NS-NS 
binaries and ^OMpc for a BH-NS binary with a mass ratio of 4 : 1, and a dedicated merger 
search on this data did not find any evidence for GW candidates [2]. 

2.4 Bursting magnetars 

Soft-gamma repeaters (SGRs) and anomalous X-ray pulsars are X-ray pulsars with quiescent 
soft (2 — lOkeV) periodic X-ray emissions with periods ranging from 5 to 10 s. They exhibit 
repetitive outbursts lasting ~0.1s which reach peak luminosities of ^10 42 ergs' 1 in X-rays 
and y-rays. There are a number of known SGRs and anomalous X-ray pulsars [86, 106], some 
of which have had rare hard spectrum giant flares with luminosities of up to 10 47 erg s _1 . The 
favoured model for these objects is a magnetar, a neutron star with an extreme magnetic field 
oiB ~ 10 15 G. Giant flares are believed to be caused either by magnetic stresses fracturing the 
magnetar crust and leading to a large-scale rearrangement of the internal field [139] or by a 
large-scale rearrangement of the magnetospheric field due to magnetic reconnection [70, 102]. 
The sudden release of energy and magnetic field rearrangement lead to the creation and 
acceleration of pair plasma that may have some baryon loading, thus leading to the emission 
of HENs in py reactions [78]. The detectability of the 2004 giant are of the Galactic SGR 
1806-20 by HEN detectors was estimated in [88]: detectors such as IceCube and ANTARES 
should detect multiple HEN events from similar Galactic SGR eruptions, provided the baryon 
loading is sufficiently high. The AMANDA-II detector, which was operating during the 
giant flare of SGR 1806-20, did not detect HENs [23]. A search of IceCube data for HENs 
from regular (non-giant) Galactic SGR flares also found no significant coincident events [13]. 
Estimates based on [88] for the 5-line ANTARES detector show that, d^o ss 200 kpc for 
baryon-rich flares, suggesting that similar flares could be detected from anywhere within the 
Galaxy. 

Significant emission of GWs in SGR giant flares may come from the potential excitation 
of nonradial pulsational modes with kHz- frequencies in the magnetar [62] . Theoretical upper 
limits on the possible strength of the GW emission were placed by [87] and [56], based on the 
the energy reservoir associated with a change in the magnetic potential energy of the mag- 
netar. They found an upper limit for the emitted GW energy of 10“' Mqc 2 to 10“ 6 M 0 c 2 , 
which can be probed by the LIGO/Virgo network for a Galactic source [4, 15]. However, 
studies that investigated the excitation of magnetar pulsational modes in more detail sug- 
gest much weaker emission that may not be detectable even with advanced-generation GW 
observatories [98, 153]. 

2.5 Cosmic string kinks and cusps 

Cosmic strings are topological defects that may form in the early Universe and are pre- 
dicted by grand unified theories and superstring theory [e.g., 91, 125]. Cosmic strings form 
initially as smooth loops, but through interactions and self-interactions may develop kinks 
and cusps [e.g., 125]. The kinks and cusps decay, which is expected to lead to ultra-high 
energy cosmic ray emission with energies in excess of ~ 10 11 GeV and up to the Planck 
scale [47, 82], including ultra-high-energy neutrinos (UHENs; e.g., [37, 45, 101]) and GW 
bursts [e.g., 59, 60, 111, 134]. 

While not designed specifically for UHENs, HEN detectors like ANTARES and IceCube 
have some sensitivity to UHENs in the > 10 11 GeV energy range. Up to a few events per 
year for a km ! -scale detector are predicted in [35], depending on details of the underlying 
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emission model. Since Earth is opaque to UHENs, downgoing events must be selected. Since 
we are only considering ANTARES data for neutrinos that have passed through Earth (see 
section 3.1), the present data set does not contain any potential UHEN events. 

The search for very energetic HENs performed with one year of IceCube-40 data did not 
reveal any neutrinos in the 10 6 GeV to 10 9 GeV energy range [9], but no limits on UHENs 
were reported. A number of dedicated UHEN experiments exist, including ANITA [71], 
NuMoon [132] and others, but have not yet constrained many emission scenarios from cosmic 
strings [e.g., 101]. 

The rate of GW bursts from a network of cosmic strings depends on the string tension 
and other network parameters, and individual bursts may be detectable with advanced de- 
tectors [60, 134]. The burst shape is expected to be generic, so that matched-filtering GW 
analysis approaches may be employed. A first search for GW bursts from cosmic string cusps 
in 15 days of LIGO data from early 2005 did not reveal any candidate events [17]. 

3 GW and HEN detectors 

3.1 The ANTARES neutrino telescope 

Since the Earth acts as a shield against all particles except neutrinos, a neutrino telescope 
mainly uses the detection of upgoing muons as a signature of muon-neutrino charged-current 
interactions in the matter around the detector. The ANTARES detector (Astronomy with 
a Neutrino Telescope and Abyss environmental RESearch) is currently the only deep sea 
higli-energy-neutrino telescope and is operating in the Northern hemisphere [28]. The tele- 
scope covers an area of about 0.1km 2 on the sea bed, at a depth of 2475 m, 40 km off the 
coast of Toulon, France. The detector is a three-dimensional array of photomultiplier tubes 
(PMTs) [30], hosted in pressure resistant glass spheres, called optical modules (OMs) [36]. In 
its full configuration, it is composed of 12 detection lines, each comprising up to 25 triplets of 
PMTs, storeys, regularly distributed along 350 m, the first storey being located 100 m above 
the sea bed. The first detection line was installed and connected in early 2006; the second 
line was put in operation in September 2006 and three more lines were connected in January 

2007, so that a total of 5 lines were taking data in 2007. Five additional lines, together 
with an instrumentation line (containing an ensemble of oceanographic sensors dedicated to 
the measurement of environmental parameters), were connected by the end of 2007. The 
telescope reached its nominal configuration, with 12 lines immersed and taking data, in May 

2008. 

The three-dimensional grid of PMTs is used to measure the arrival time and position of 
Cherenkov photons induced by the passage of relativistic charged particles through the sea 
water. This information, together with the characteristic emission angle of the light (about 
43 degrees), is used to determine the direction of the muon and hence infer that of the 
incident neutrino. The accuracy of the direction information allows to distinguish upgoing 
muons, produced by neutrinos, from the overwhelming background from downgoing muons, 
produced by cosmic ray interactions in the atmosphere above the detector [32]. Installing 
the detector at great depths serves to attenuate this background and also allows to operate 
the PMTs in a dark environment. At high energies the large muon range makes the sensitive 
volume of the detector significantly greater than the instrumented volume. By searching for 
upgoing muons, the total ANTARES sky coverage is 3.57T sr, with most of the Galactic plane 
being observable and the Galactic Center being visible 70% of the sidereal day. 
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3.2 Network of interferometers 


3.2.1 LIGO 

LIGO is a network of interferometric gravitational wave detectors consisting of three inter- 
ferometers in the U.S.A.. These detectors are all kilometer-scale power-recycled Michelson 
laser interferometers with orthogonal Fabry-Perot arms [18] able to detect the quadrupolar 
strain in space produced by the GW. Multiple reflections between mirrors located at the end 
points of each arm extend the effective optical length of each arm, and enhance the sensitivity 
of the instrument. 

There are two LIGO observatories: one located at Hanford, WA and the other at 
Livingston, LA. The Hanford site houses two interferometers: one with 4 km arms, denoted 
HI, and a second with 2km arms, denoted H2. The Livingston observatory has one 4km 
interferometer, LI. The observatories are separated by a distance of 3000 km, corresponding 
to a time-of-flight separation of 10 ms. 

The LIGO instruments are designed to detect gravitational waves with frequencies 
ranging from ~ 40 Hz to several kHz, with a maximum sensitivity near 150 Hz (see figure 1). 
In fact, seismic noise dominates at lower frequencies and the sensitivity at intermediate 
frequencies is determined mainly by thermal noise, with contributions from other sources. 
Above ~ 200 Hz, laser shot noise corrected for the Fabry-Perot cavity response yields an 
effective strain noise that rises linearly with frequency. The average sensitivities of the HI 
and LI detectors during the second year of the S5 run were about 20% better than the 
first-year averages, while the H2 detector had about the same average sensitivity in both 
years. 


3.2.2 Virgo 

The Virgo detector, VI, is in Cascina near Pisa, Italy. It is a 3 km long power- recycled Michel- 
son interferometer with orthogonal Fabry Perot arms [21]. The main instrumental difference 
with respect to LIGO is the seismic isolation system based on super-attenuators [48], chains 
of passive attenuators capable of filtering seismic disturbances. The benefit from super- 
attenuators is a significant reduction of the detector noise at very low frequency (<40Hz) 
where Virgo surpasses the LIGO sensitivity. During 2007, above 300 Hz, the Virgo detector 
had sensitivity similar to the LIGO 4 km interferometers, while above 500 Hz it is dominated 
by shot noise, see figure 1. 

The time-of-flight separation between the Virgo and Hanford observatories is 27 ms, 
and 25 ms between Virgo and Livingston. Due to the different orientation of its arms, the 
angular sensitivity of Virgo is complementary to that of the LIGO detectors, Virgo therefore 
enhances the sky coverage of the network. Moreover, simultaneous observations of multiple 
detectors are crucial to reject environmental and instrumental effects. 

At the time of writing the LIGO and Virgo interferometers are undergoing upgrades to 
“advanced” configurations with distance sensitivity improved by approximately a factor of 
10 [79]. The advanced instruments will commence operations around 2015. 

3.3 Joint data taking periods 

The fifth LIGO science run, S5 [14], was held from 2005 November 4 to 2007 October 1. 
Over one year of science-quality data were collected with all three LIGO interferometers in 
simultaneous operation at their design sensitivity, with duty factors of 75%, 76%, and 65% 
for HI, H2, and LI. The Virgo detector started its first science run, VSR1 [22], on 2007 
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Figure 1. Noise amplitude spectral densities of the four LIGO and Virgo detectors during S5. 


May 18. The Virgo duty factor over VSR1 was 78%. During this period, ANTARES was 
operating in its 5 line configuration. The concomitant set of ANTARES 5- line (5L), VSR1 
and S5 data covers the period between January 27 and September 30, 2007; these data are 
the subject of the analysis presented here. 

4 Selection of HEN candidates 

4.1 HEN data sample 

The ANTARES data sample used in the analysis is composed of runs from 2007 selected 
according to various quality criteria, based mainly on environmental parameters (e.g. sea cur- 
rent, counting rates), configuration and behaviour of the detector during the given run (e.g. 
duration of the run, alignment of the detector). Two basic quantities are used to characterise 
the counting rate of a given OM: the baseline rate ( 40 K activity and bioluminescence) and 
the burst fraction (flashes of light emitted by marine organisms). The baseline rate repre- 
sents the most probable counting rate of a given OM computed from the rate distributions in 
each PMT over the whole run (typically a few hours). The burst fraction corresponds to the 
fraction of time during which the OM counting rates exceed by more than 20% the estimated 
baseline. The data selected for this search are required to have a baseline rate below 120 kHz 
and a burst fraction lower than 40%, with 80% of all OMs being active. With these quality 
criteria, the active time is 103.4 days out of the 244.8 days of the 5-line period. Finally, when 
restricting the data to the concomitant period with LIGO/Virgo, the remaining equivalent 
time of observation is T 0 b s = 91 days. 

4.2 Trigger levels 

The ANTARES trigger system is multi-level [31]. The first level is applied in situ, while the 
remaining levels intervene after all data are sent to the shore station and before being written 
on disk. Trigger decisions are based on calculations done at three levels. The first trigger 
level, L0, is a simple threshold of about 0.3 photo-electron (pe) equivalent charge applied to 
the analog signal of the PMT. The second level trigger, LI, is based on two coincident L0 hits 
in the same storey within 20 ns or hits with large charge (> 3pe or 10 pe depending on the 
configuration). The L2 trigger requires the presence of at least five LI hits in a 2.2 /rs time 
window (roughly the maximum muon transit time across the detector) and that each pair of 
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LI hits are causally related according to the following condition: A tij < dy7i/c + 20ns. Here 
A tij and are the time difference and distance between hits i and j, c is the speed of light 
in vacuum and n is the index of refraction. 


4.3 Reconstruction strategy 

Hits selected according to the criteria described in section 4.2 are then combined to recon- 
struct tracks using their arrival time and charge as measured by the corresponding OM. 
Muons are assumed to cross the detector at the speed of light along a straight line from 
which the induced Cherenkov light originates. The time and charge information of the hits 
in the PMTs is used in a minimisation procedure to obtain the track parameters, namely, 
its direction and the position (xo,yo,zo) of one track point at a given time to- The 

reconstruction algorithm used for this analysis is a fast and robust method [34] which was 
primarily designed to be used on-line. 


4.3.1 Description of the algorithm 

The algorithm is based on a ^-minimisation approach. Its strict hit selection leads to a 
high purity up-down separation while keeping a good efficiency. The exact geometry of the 
detector is ignored: the detector lines are treated as straight and the 3 OMs of each storey 
are considered as a single OM centered on the line. Thus, the hit’s altitude corresponds 
to the optical modules altitude. All hits at the same floor in coincidence within 20 ns are 
merged into one hit. The time of the merged hit is that of the earliest hit in the group and 
its charge is the sum of the charges. The algorithm uses the LI hits as a seed for the hit 
selection. It requests a coincidence of 2 LI hits in two adjacent floors within 80 ns or 160 ns 
in two next-to- adjacent floors. The quality of the reconstruction is measured by a y 2 -like 
variable with NDF degrees of freedom, based on the time differences between the hit times t, 
and the expected arrival time tj of photons from the track or bright-point (see section 4.4). 
The quality function is then extended with a term that accounts for measured hit charges qi 
and the calculated photon travel distances c (- : 


1 y 2 ' At 2 Q(qi) D jd]) 

NDF ^ of q d 0 
1=1 L L - 1 


(4.1) 


In this expression, cq is the timing error, set to 10 ns for charges larger than 2.5 pe and to 
20 ns otherwise. Af * = tj — t, is the time residuals between the hit time t t and the expected 
arrival time of the photons tj from the muon track. In the second term, q is the average 
hit charge calculated from all hits which have been selected for the fit and do — 50 m is the 
typical distance at which the signal in one PMT from a Cherenkov light front is of the order 
of 1 pe. The function Q(qi) accounts for the angular acceptance of the OMs, while D(dJ) 
penalises large amplitude hits originating from large distance tracks. A proper cut on the fit 
quality parameter allows the isolation of a high purity neutrino sample, which is crucial in 
the subsequent analysis. 


4.3.2 Azimuthal degeneracy of the reconstruction 

For a particle trajectory reconstructed from a Cherenkov cone giving hits on only two straight 
detector lines, there always exists an alternative trajectory having an identical y 2 value, but 
a different direction. The degenerate trajectory is the mirror image of the original track in 
the plane formed by the two lines. As a consequence, each event reconstructed with only two 
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Figure 2. Distribution of the sine of the declination 5 of selected events (black points), compared to 
Monte-Carlo expectations (sum of atmospheric muons and atmospheric neutrinos, orange (or grey) 
points). 


lines will have two equiprobable arrival directions, which must be taken into account during 
the subsequent GW analysis. 

4.4 Criteria for HEN event selection 

The initial sample of reconstructed events contains both upgoing neutrino induced muons and 
downgoing muons from cosmic ray interactions in the atmosphere. Some of the atmospheric 
muons are misreconstructed as upgoing and the selection cuts, based on Monte-Carlo simu- 
lations, are devised to reduce this contamination so as to maximise the discovery potential. 
A minimum of 6 hits on at least 2 lines are required to reconstruct a track. Only upgoing 
tracks are kept for further analysis. Quality cuts are then applied based on two quantities 
computed according to equation (4.1). The first parameter used, Xt> quality factor 

associated with the reconstructed particle track, whereas the second one, Xb> is associated 
with a bright-point, light emitted from a point-like source inside the detector. This rejects 
events from large electromagnetic showers, likely to appear in downgoing muon bundles for 
instance. 

A cut on xt reduces the number of such events and decreases the contribution of mis- 
reconstructed muons in the background. Further cuts are applied on Xt depending on the 
arrival direction of the candidate - the muon contamination increases close to the horizon - 
which reduce the fraction of misreconstructed muons to less than 20% over the whole sample, 
while optimising the sensitivity (see section 4.6 and [77]). 

Figure 2 shows the distribution of the sine of the declination of the events selected with 
the final cuts, which is globally consistent with background. 

4.5 Angular error 

The distribution of the space angle fl between the true neutrino arrival direction and the 
reconstructed muon track can be described by a log-normal distribution: 



\/2n (O — #o)oo 

where 6q is a location parameter, op is related to the shape of the distribution and m-o is 
a scaling parameter. In all cases for our study, the location parameter 8q is close to zero, 


(4.2) 
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Figure 3. Example of space angle distribution with the associated fit to equation (4.2) obtained with 
a sample of Monte- Carlo HEN events, for a given declination and a given number of hits. The arrows 
indicate the 50 th (median) and the 90 th percentiles of the distribution. The distribution is normalised 
to unity. 


and (fi — #o) > 0 is always satisfied. This distribution depends on the energy associated 
to the track (estimated through the number of photons detected) and its declination. This 
parametrisation is used during the GW search to compute the significance of a hypotheti- 
cal signal for the scanned directions inside the angular search window centred around the 
reconstructed neutrino arrival direction. Figure 3 shows an example of distribution of the 
space angle for a sample of Monte Carlo neutrinos with an E~ 2 spectrum, together with the 
best-fit parametrisation and the 50 th and 90 th percentiles of the distribution. 

One of the main variables to describe the performance of a neutrino telescope is the 
angular resolution, defined as the median of the distribution of the angle between the true 
neutrino direction and the reconstructed track, also indicated in figure 3. This number is 
estimated from simulations. 

For those events of our selected sample reconstructed with at least three lines the angular 
resolution is, assuming an E~ 2 energy spectrum, ~ 2.5° at 100 GeV, improving to 1° around 
100 TeV. For 2-line events, when selecting the reconstructed track closer to the true direction, 
the angular accuracy varies between 3° at low energy (100 GeV) and 2.5° at high energy 
(100 TeV). 

We define the angular search window for the GW analysis as the 90 th percentile of the 
distribution, also indicated in figure 3; this window lies between 5° and 10° for 3-line events, 
depending on declination, and between 10° and 15° for 2-line events. 

We note that the typical angular distance between galaxies within lOMpc is a few 
degrees [148], much smaller than the typical size of the 90 th percentile error region for our 
HEN events. This implies that we can associate a potential host galaxy to any of the HEN 
candidates if it turns out to be of cosmic origin. 

4.6 Analysis sensitivity and selected HEN candidates 

The limit-setting potential of the analysis, or sensitivity, has been quantified for the whole 
5 line data period. Specifically, the sensitivity is defined as the median 90% upper limit 
obtained over an ensemble of simulated experiments with no true signal. The sensitivity 
depends on the declination of the potential source. For our sample and assuming an E~ 2 
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7 


Sky Map 



Figure 4. Skymap of the selected 216 HEN events in equatorial coordinates. A line connects the 
associated mirror solutions for events reconstructed with two lines as described in section 4.3.2. 


steady flux, using the selection criteria described, the best sensitivity has been estimated 
to be E ~ 10 -6 GeV cm -2 s -1 . This best sensitivity is reached below —47°; i.e., at 
declinations which are always below the horizon at the latitude of ANTARES (43°N). 

With the selection previously described, 181 runs corresponding to 104 days of live time 
were kept for the analysis. The selection has been divided into events reconstructed with 

2 lines and events with at least 3 lines. Each of the mirror solutions for 2 line events will 
be searched for possible counterparts in the subsequent GW analysis. This results in 216 
neutrinos to be analysed: 198 with two possible directions and 18 reconstructed with at least 

3 lines. Figure 4 is a sky map of the candidate HEN events, where the degenerate solutions 
for 2 line events can be seen. 

Of these HEN events, 158 occurred at times when at least two gravitational-wave de- 
tectors were operating. Since two or more detectors are required to discriminate GW signals 
from background noise (as described in section 5.2), in the following we consider only these 
remaining 158 HEN candidates: 144 2- line events and 14 3-line events. 1 

Finally, we note that IceCube operated in its 22-string configuration for part of 2007 [13]. 
However, this data was only used for time-dependent searches applied to source directions 
with observed X-ray or gamma-ray emission, such as GRBs; there were no untriggered, time- 
dependent searches over the sky. Furthermore, a comparison of ANTARES and IceCube 
sensitivities in 2007 indicates that the bulk of our HEN neutrino triggers come from declina- 
tions (the southern sky) such that it is unlikely that IceCube could have detected the source 
independently. 


1 Details of each of the HEN candidate events are given at https://dcc.ligo.org/cgi-bin/DocDB/ShowDocu 
ment?docid=pl 200006 . 
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5 GW search method 


5.1 Search procedure 

One of the simplest searches that may be performed combining GW and HEN data is a 
triggered analysis that scans GW data around the time of the putative neutrino event by 
cross-correlating data from pairs of detectors. This search exploits knowledge of the time 
and direction of the neutrino event to improve the GW search sensitivity. We use the X- 
Pipeline algorithm [137], which has been used in similar searches for GWs associated with 
GRBs [6, 20]. X-PlPELlNE performs a coherent analysis of data from arbitrary networks 
of gravitational wave detectors, while being robust against background noise fluctuations. 
Each trigger is analysed independently of the others, with the analysis parameters optimised 
based on background noise characteristics and detector performance at the time of that 
trigger, thereby maximising the search sensitivity. 

5.2 GW event analysis 

In our GW search, a neutrino candidate event is characterised by its arrival time, direction, 
and angular search window (and mirror- image window, for the 2-line events). Also important 
is the range of possible time delays (both positive and negative) between the neutrino emission 
and the associated gravitational-wave emission. This quantity is referred to as the on-source 
window for the neutrino; it is the time interval which is searched for GW signals. We use a 
symmetric on-source window of ±496 s [41] , which is conservative enough to encompass most 
theoretical models of GW and HEN emission. The maximum expected time delay between 
GWs and HENs due to a non-zero mass effect for either particle is much smaller than the 
coincidence windows used. 

The basic search procedure follows that used in [20] . All detectors operating at the time 
of the trigger and which pass data-quality requirements are used for the GW search. The 
data from each detector are first whitened and time-delayed according to the sky location 
being analysed so that a GW signal from that direction would appear simultaneous in each 
data stream. The data are then Fourier transformed to produce time-frequency maps. The 
maps are summed coherently (using amplitude and phase) with weighting determined by 
each detector’s frequency-dependent sensitivity and response to the sky location in question; 
the weightings are chosen to maximise the signal-to-noise ratio expected for a circularly 
polarized GW signal, 2 which is the expected polarisation for a GW source observed from 
near the rotational axis [94]. A threshold is placed on the map to retain the largest 1% 
of pixels by energy (squared amplitude). Surviving pixels are grouped using next-nearest- 
neighbours clustering; each cluster of pixels is considered as a candidate GW event. The 
event cluster is assigned a combined energy by summing the energy values of its constituent 
pixels; this combined energy is used as the ranking statistic for the events. 

In addition to the marginalised circular polarization sum, a second ranking statistic 
is computed based on a maximum-likelihood analysis of the event assuming power-law dis- 
tributed background noise with no assumption on the GW polarization. In practice this 
statistic is often found to provide signal-noise separation due to the non-Gaussian nature of 
the GW detector noise. Other combinations of the data are also constructed. Of particular 
importance are “null” combinations designed to cancel out the GW signal from the given 

Empirically it is found that the circular polarisation restriction also improves the overall detection prob- 
ability for linearly polarised GWs, as the resulting background reduction outweighs the impact of rejecting 
some linearly polarised GWs. 
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sky location; comparison to corresponding “incoherent” combinations provides powerful tests 
for identifying events due to background noise fluctuations [53], and are described in detail 
in [145]. Events are also characterised by their duration, central time, bandwidth, and central 
frequency. 

The time-frequency analysis is repeated for Fourier transform lengths of 1/128, 1/64, 
1/32, 1/16, 1/8, 1/4 s, to maximise the sensitivity to GW signals of different durations. It is 
also repeated over a grid of sky positions covering the 90% containment region of the HEN. 
This grid is designed such that the maximum relative timing error between any pair of GW 
detectors is less that 0.5 ms. When GW events from different Fourier transforms lengths 
or sky positions overlap in time-frequency, the highest-ranked event is kept and the others 
discarded. Finally, the events are decimated to a rate of 0.25 Hz before being written to disk. 

This time-frequency analysis is performed for all of the data in the ± 496 s on-source 
window. To estimate the significance of the resulting GW candidates, the same analysis is 
repeated for all coincident data in the off- source window, defined as all data within ±1.5 hours 
of the neutrino time, excluding the on-source interval. The same set of detectors and data- 
quality requirements as in the on-source analysis are used for the off-source data. These off- 
source data provide a sample of background that does not contain any signal associated with 
the neutrino event, but with statistical features similar to the data searched in association 
with the neutrino. To enlarge the background sample, we also repeat the off-source analysis 
after applying time shifts of multiples of 6 s to the data from one or more detectors; with 
such time slides we were able to produce O(10 3 ) background trials for each HEN. 

Finally, the analysis is repeated after “injecting” (adding) simulated GW signals to the 
on-source data. The amplitudes and morphologies tested are discussed in section 6.3.1. We 
use these simulations to optimise and assess the sensitivity of the search, as discussed below. 


5.3 GW search optimisation 

The sensitivity of searches for gravitational-wave bursts tends to be limited by the presence 
of non-Gaussian fluctuations of the background noise, known as glitches. To reduce this 
background, events that overlap in time within known instrumental and/or environmental 
disturbances are discarded. In addition to this “veto” step, GW consistency tests comparing 
the coherent and incoherent energies are applied to each event [145]. These tests are applied 
to the on-source, off-source and injection events; events failing one or more of these tests are 
discarded. The thresholds are optimised by testing a preset range of thresholds and selecting 
those which give the best overall detection efficiency at a fixed false alarm probability of 
1% when applied to a random sample of background and injection events (the on-source 
events are not used; i.e. , this is a blind analysis). These tests also determine which of the 
two ranking statistics discussed in section 5.2 (based on circularly polarized GW energy 
or powerlaw noise) gives the better detection efficiency; the winner is selected as the final 
ranking statistic. 

Once the thresholds have been fixed, these consistency tests are applied to the on-source 
events and to the remaining off-source and injection events (those not used for tuning). The 
surviving on-source event with the largest significance (highest energy or powerlaw statistic) 
is taken to be the best candidate for a gravitational wave signal and is referred to as the 
loudest event [49]. All surviving on-source events are assigned a false alarm probability by 
comparison to the distribution of loudest events from the off-source trials. Any on-source 
event with probability p < 0.01 is subjected to additional checks to try to determine the 
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origin of the event and additional background time slide trials are performed to improve the 
accuracy of the false alarm probability estimate. 

After the p values have been determined for the loudest events associated with each 
of the 158 HEN events, the collective set of p values is tested for consistency with the null 
hypothesis (no GW signal) using the binomial test, discussed in section 6.2. We also set a 
frequentist upper limit on the strength of gravitational waves associated with each neutrino 
trigger, as discussed in section 6.3. 

5.4 Low-frequency and high-frequency GW analyses 

Given our knowledge of possible GW sources discussed in section 2, the most likely detectable 
signals at extra-galactic distances are in the low-frequency band (/ < 500 Hz), where our 
detectors have maximum sensitivity, see figure 1. At the same time, the computational cost 
of the X-Pipeline analysis increases at high frequencies. This is due in part to the extra data 
to be analysed, but also to the need for finer-resolution sky grids to keep time delay errors 
much smaller than one GW period. We therefore split the gravitational wave band into two 
regions: 60 Hz to 500 Hz and 500 Hz to 2000 Hz. The low-frequency band is analysed for all 
HEN events — such a search is computationally feasible while covering the highest-sensitivity 
region of the GW detectors. However, compact objects such as neutron stars or collapsar 
cores have characteristic frequencies for GW emission above 500 Hz. Such emissions might 
be detectable from Galactic sources such as soft gamma repeater giant flares, or possibly 
from nearby galaxies. Since the computational cost of a high-frequency search for all HEN 
events is prohibitive with the current analysis pipeline, we perform the 500 Hz to 2000 Hz 
analysis on the 3- line HEN events only. The 3-line events are a small subset (^10%) of the 
total trigger list and have the smallest sky position uncertainties, and therefore the smallest 
computational cost for processing. To reduce the computational cost further, we use the same 
sky grid for the high-frequency search as was used at low frequencies, after determining that 
the loss of sensitivity is acceptable. The high-frequency analysis is performed independently 
of the low-frequency analysis (independent tuning, background estimation, etc.) using the 
identical automated procedure. In the following sections we will present the results of the 
low- frequency and high-frequency searches separately. 

6 Coincident search results 
6.1 Per-HEN GW candidates 

We analysed GW data in coincidence with 158 neutrino candidates for the low frequency 
search, and 14 neutrino events for the high frequency search. In the low frequency analy- 
sis, only one neutrino trigger had a corresponding GW event with false alarm probability 
below the threshold of p — 0.01 to become a candidate event. We found no candidates in 
the high frequency search. For the low-frequency candidate, additional time shifts totaling 
18064 background trials yielded a refined false alarm probability of p — 0.004, which is not 
significant given a trials factor of 158 (this statement is quantified below). This event came 
from analysis of the HI, H2, and VI data; follow-up checks were performed, including checks 
of detector performance at the time as indicated by monitoring programs and operator logs, 
and scans of data from detector and environmental monitoring equipment to look for anoma- 
lous behaviour. While these checks did not uncover a physical cause for the event, they did 
reveal that it occurred during a glitching period in VI. We conclude that we have no clear 
gravitational wave burst signal associated with any of our sample of 158 neutrino events. 
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Figure 5. Distribution of observed p values for the loudest GW event associated with each neutrino 
analysed in the low frequency analysis. The red dot indicates the largest deviation of the low p tail 
from the uniform distribution null hypothesis; this occurs due to having the three loudest events below 
P 3 ~ 0.013. Deviations this large or larger occur in approximately 64% of experiments under the null 
hypothesis. The black line shows the threshold for a 5-sigma deviation from the null hypothesis. 


6.2 Search for a cumulative excess: binomial test 

A quantitative analysis of the significance of any candidate gravitational-wave event must 
take account of the trials factor due to the number of neutrino events analysed. We use 
the binomial test, which has been applied in previous GRB-triggered GW searches [16, 20]. 
Under the null hypothesis, the false alarm probabilities p for each HEN loudest event are 
expected to be uniformly distributed between 0 and 1. The binomial test compares the 
measured p values to the null distribution to determine if there is a statistically significant 
excess of (one or more) small p values which may be due to gravitational wave signals. 

Briefly, the binomial test sorts the set of N measured loudest event probabilities in 
ascending order: pi < P2 < P3 < ■ ■ ■ < Pn- For each i E [1, Wail] we compute the binomial 
probability P>i{pi ) of getting i or more events with p values < pp 

N i\j\ 

- E (at -tjiti rfa -«)"-*■ i 6 - 1 ) 

Here N is the number of HEN events analysed (158 in the 60 Hz to 500 Hz band and 14 in 
the 500 Hz to 2000 Hz band), and -/V ta ii is the number of the smallest p values we wish to test. 
We choose Wail to be 5% of AT; i.e., Wail — 8 for the low frequency band and Wail = 1 for 
the high frequency band. 

The lowest P>i(pi) for i E [1, Wail] is taken as the most significant deviation from the 
null hypothesis. To assess the significance of the deviation, we repeat the test using p values 
drawn from a uniform distribution and count the fraction of such trials which give a lowest 
P>i(jPi) smaller than that computed from the true measured p values. 

Figures 5 and 6 show the cumulative distribution of p values measured in the low- and 
high-frequency analyses. In both cases the measured p values are consistent with the null 
hypothesis. 
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Figure 6. Distribution of observed p values for the loudest GW event associated with each neutrino 
analysed in the high frequency analysis. The red dot indicates the largest deviation of the low p tail 
from the uniform distribution null hypothesis; since iV ta ii = 1, this is constrained to occur for p\. 
Deviations this large or larger occur in approximately 66% of experiments under the null hypothesis. 
The black dot shows the threshold for a 5-sigma deviation from the null hypothesis. 


6.3 GW upper limits 

The sensitivity of the GW search is determined by a Monte-Carlo analysis. For each neutrino 
trigger, we add simulated GW signals to the on-source data and repeat the analysis described 
in section 5.2. We consider a simulated signal detected if it produces an event louder than the 
loudest on-source event after all event tests have been applied. We define a 90% confidence 
level lower limit on the distance to the source as the maximum distance Dg 0 % such that for 
any distance D < Dg 0 % the probability of detection is 0.9 or greater. 


6.3.1 Injected waveforms 

As in GRB-triggered searches, we use a mix of ad hoc and astrophysically motivated GW 
waveforms. The ad hoc waveforms are Gaussian-modulated sinusoids: 
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Here /o is the central frequency, to is the central time, and r is the duration parameter. 
This waveform is consistent with the GW emission from a rotating system viewed from an 
inclination angle l to the rotational axis. We select the inclination uniformly in cos t with 
l 6 [0°,5°]. This corresponds to a nearly on-axis system, such as would be expected for 
association with an observed long GRB. We chose r — l//o, and use central frequencies of 
100 Hz, 150 Hz, and 300 Hz for the low-frequency analysis and 554 Hz and 1000 Hz for the 
high-frequency search. The quantity h ISS is the root-sum-square signal amplitude: 


h 


rss — 



{h\(t) + h x (f)) dt . 


(6.4) 
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For the small values of l considered here (t < 5°) this amplitude is related to the total energy 
i?GW in a narrow-band gravitational-wave burst by 

fi GW — ^—Q- h rssfoD 2 . (6.5) 

For astrophysical injections we use the gravitational-wave emission of inspiraling neutron 
star and black hole binaries, which are widely thought to be the progenitors of short GRBs. 
Specifically, we use the post-Newtonian model for the inspiral of a double neutron star system 
with component masses mi — m2 — 1.35 Mq, and the one for a black-hole - neutron-star 
system with m\ — 5 Mg, m 2 — 1.35Mg. We set the component spins to zero in each case. 
Motivated by estimates of the jet opening angle for short GRBs, we select the inclination 
uniformly in cost with t 6 [0°,30°]. 

For each HEN trigger, the injections are distributed uniformly in time over the on-source 
window. The injection sky positions are selected randomly following the estimated proba- 
bility distribution (4.2) for the HEN trigger, to account for the uncertainty in the true HEN 
direction of incidence. The polarization angle (orientation of the rotational axis on the sky) is 
distributed uniformly. Finally, the amplitude and arrival time at each detector is perturbed 
randomly to simulate the effect of calibration errors in the LIGO and Virgo detectors. 

6.3.2 Exclusion distances 

For each waveform type we set a 90% confidence level lower limit on the distance to a GW 
source associated with a given HEN trigger. 3 This is defined as the maximum distance Dg 0 % 
such that for any distance D < D m % there is a probability of at least 0.9 that such a GW sig- 
nal would have produced an event louder than the loudest on-source event actually measured. 
For inspirals, each distance corresponds to a well-defined amplitude. We can associate an 
amplitude to each distance for the sine-Gaussian waveforms as well, by assuming a fixed en- 
ergy in gravitational waves. For concreteness, we select £gw — 10 ~ 2 Mqc 2 . This corresponds 
to the optimistic limit of possible gravitational- wave emission by various processes in the col- 
lapsing cores of rapidly rotating massive stars ([66, 67, 93, 124], and discussion in section 2); 
more conservative estimates based on simulations have been made in [64, 118, 119, 131, 138]. 

For each type of gravitational wave simulated, the distributions of exclusion distances 
for our neutrino sample are shown in figures 7 and 8. For binary neutron star systems of 
(1.35 — 1.35)Mq and black hole - neutron star systems of (5 — 1.35 )M 0 typical distance 
limits are 5 Mpc and 10 Mpc respectively. For the sine-Gaussian waveforms with £gw — 
10“ 2 Mqc 2 we find typical distance limits between 5 Mpc and 17 Mpc in the low-frequency 
band and of order 1 Mpc in the high-frequency band. For other Eqw the limits scale as 
Dq 0% oc (£'gw/10“ 2 Mqc 2 ) 1 / 2 . For example, for £gw = 10~ h Mqc 2 (typical of core-collapse 
supernovae) a signal would only be observable from a Galactic source. 

7 Astrophysical implications 

Observational constraints on joint sources of GW and HEN signals have been derived in [43]. 
However, they are based on the interpretation and the combination of previously published 
and independent GW and HEN observational results. The results presented in this section 
are the first derived from a joint GW-HEN analysis, using concomitant data obtained with 
LIGO/Virgo and ANTARES. 

3 Upper limits for each waveform and HEN trigger are available at https://dcc.ligo.org/cgi-bin/DocDB/ 
ShowDocument?docid=pl200006. 
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Figure 7. Low-frequency analysis: the top plot is the histogram for the sample of analysed neutrinos 
of the distance exclusions at the 90% confidence level for the 3 types of sine-Gaussian models consid- 
ered: 100 Hz, 150 Hz and 300 Hz. A standard siren gravitational wave emission of Egw — 10 -2 M 0 c 2 
is assumed. The bottom plot shows the distance exclusions for the 2 families of binary inspiral models 
considered: NS-NS and BH-NS. 



Figure 8. High-frequency analysis: the histogram for the sample of analysed neutrinos of the distance 
exclusions at the 90% confidence level for the 2 frequencies of circular sine-Gaussian models considered: 
554 Hz and 1000 Hz. 


7.1 Upper limits on GW-HEN populations 

The present search for GW and HEN correlations in space and time revealed no evidence for 
coincident events. This implies a 90% confidence level upper limit on the rate of detectable 
coincidences of 2.3/T 0 b s , where T 0 b s ~ 90 days is the duration of coincident observations. 
This can be expressed as a limit on the rate density (number per unit time per unit volume) 
PGW-HEN of joint GW-HEN sources: 


PGW-HEN < 


2.3J h 
VT ohs ‘ 


(7.1) 
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Here J-b is the beaming factor (the ratio of the total number of sources to the number with 
jets oriented towards Earth 4 ), and V is the volume of universe probed by the present analysis 
for typical GW-HEN sources. 

We take as fiducial sources two classes of objects: the final merger phase of the coa- 
lescence of two compact objects (short GRB-like), or the collapse of a massive object (long 
GRB-like), both followed by the emission of a relativistic hadronic jet. We define the HEN 
horizon as the distance for which the probability to detect at least 1 HEN in ANTARES with 
5 lines is 50%. In the case of short GRBs (SGRBs), the HEN horizon is estimated to be 
d§o — 4Mpc using [76], while the typical GW horizon from the inspiral model is 5Mpc 
to lOMpc depending on the binary masses. For long GRBs (LGRB) the HEN horizon 
increases to d^o — 12Mpc using [76]. The GW emission associated with long GRBs is 
highly uncertain; our optimistic assumption of Eq \y = 10 _2 MqC 2 at low frequencies gives 
a typical horizon distance of 10 Mpc to 20 Mpc in GW. Using the lower of the GW and 
HEN distances in each case yields from equation (7.1) approximate limits on the popula- 
tion density. For SGRB-like sources, related to the merger of two compact objects, we find 
^GW-hen ~ Ab x 1(T 2 Mpc~ 3 yW 1 . For LGRB-like sources, related to the collapse of massive 
stars, we find Pqv^hen < ^-^aoi 2 x 10 -3 Mpc -3 yr -1 , where Eo.oi = E GW /10~ 2 M & c 2 . 

7.2 Comparison of limits with existing estimates 

After correcting for beaming effects, a local rate density of SGRBs of psgrb ~ 10 _7 Mpc _3 yr _1 
to 10 -6 Mpc _3 yr _1 is suggested in [72], [115], and [74], This is similar to the abundance of 
binary neutron star mergers, their assumed progenitors, estimated to be Pns-NS ~ 1CU 8 
MpW 3 yr~ 4 to 10“ 5 Mpc~ 3 yr~ 4 [see for example 3], and well below the reach of the present 
search (/Oq^-hen ~ x 10 -2 Mpc -3 yr -1 ). With T 0 b s — 1 yr, an improvement of a factor 10 
on the detection distance is required in order to begin constraining the fraction of mergers 
producing coincident GW— HEN signals. 

A total rate of long GRBs of plgrb ~ 3 x 10 -8 Mpc -3 yr -1 is estimated in [73], after 
correcting for beaming effects; these sources are closely related to Type II and Type Ibc 
core-collapse supernovae. The local rate of SNIbc is psNibc 2 x 10 5 Mpc 3 yr 1 [75], 
whereas psnii ~2x 10 -4 Mpc -3 yr -1 [44], relatively close to the obtained limit Pqvv-hen ~ 

— 3/2 _o _o _i 

EbE 0 0 ] x 10 1 Mpc yr under our optimistic assumptions of GW emission in this sce- 
nario. A factor 10 only is required in order to begin constraining the fraction of stellar 
collapse events producing coincident weakly beamed GW-HEN signals, which translates into 
a required improvement of 2 on the detection distance. 


8 Conclusions 

This first joint GW-HEN search using 2007 data, obtained with the ANTARES HEN tele- 
scope and the Virgo/LIGO GW interferometers, opens the way to a novel multi-messenger 
astronomy. Limits on the rate density Pgw-hen of joint GW-HEN emitting systems were 
extracted for the first time using the analysis of coincident GW-HEN data. We note that 
these limits are consistent with the ones obtained in [43] derived from independent GW- 
HEN observations. More stringent limits will be available by performing similar coincidence 
analyses using other data sets provided by the same instruments. 

4 For example, for a jet opening angle of 5° gives Tb ~ 300, while 30° gives Tb ~ 10. 
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For instance, the sixth LIGO science run S6 and second and third Virgo science 
runs VSR2,3 covered the period from 7 July 2009 to 21 October 2010. Meanwhile, the 
ANTARES telescope has taken data with first 10 then 12 active lines since the end of De- 
cember 2007. Their enhanced sensitivities should permit a combined analysis to gain the 
factor required to obtain Pgw-HEN — /°SNli/SNibc an d begin to constrain the fraction of stellar 
collapse events accompanied by the coincident emission of relativistic jets beamed towards 
Earth. The analysis of these data is underway, and a similar analysis using data from the 
LIGO/Virgo S5-VSR1 periods and the IceCube HEN telescope in its 22 string configuration 
is being finalized. 

Future observing runs involving IceCube, KM3 NeT [81], and the advanced LIGO and 
advanced Virgo projects [79], are likely to coincide as well. They will give other opportunities 
to look for potential coincident GW-HEN emissions. 
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